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We consider a model dissipative quantum-mechanical system realized by coupling a quantum 
oscillator to a semi-infinite classical string which serves as a means of energy transfer from the 
oscillator to the infinity and thus plays the role of a dissipative element. The coupling between the 
two — quantum and classical — parts of the compound system is treated in the spirit of the mean- 
field approximation and justification of the validity of such an approach is given. The equations 
of motion of the classical subsystem are solved explicitly and an effective dissipative Schrodinger 
equation for the quantum subsystem is obtained. The proposed formalism is illustrated by its 
application to two basic problems: the decay of the quasi-stationary state and the calculation of the 
nonlinear resonance line shape. 
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I. INTRODUCTION 

The rapid expansion of research frontiers into the 
nanometer and femtosecond regime demands a careful 
consideration of dissipation as a way to counterbalance 
the energy influx from the external fields in quantum me- 
chanical systems. 

The subject takes roots in the seminal work by Feyn- 
man and Vernon pj who treated a quantum object cou- 
pled to an infinite collection of oscillators as a model 
of a linear dissipative environment. By employing the 
path-integral techniques they were able to eliminate the 
environment variables and arrived at a dissipative mod- 
ification of the Green's function of the quantum object. 
This approach was further used by a number of authors, 
most notably by Caldeira and Leggett [2], who applied 
it to a specific problem of dissipative effects in tunneling 
on a macroscopic scale. A detailed account on the prob- 
lem of dissipation in quantum-mechanical systems can 
be found in a monograph by Weiss [3] while semiclassical 
approaches are also reviewed in 

Generally, the main mode of attack uses the density- 
matrix formalism 0| and the popular Lindblad technique 
[3]- This approach is relevant for open quantum systems. 
In contrast, the considerably less arduous Schrodinger- 
equation formalism is regarded to be appropriate for iso- 
lated systems that are described by a hermitian Hamil- 
tonian. 

Despite this fact, there were successful attempts to in- 
clude dissipation directly in the Schrodinger equation. 
Most often they based on the de Broglie-Bohm picture 
[3] and employed modifications to the quantum poten- 
tial [3]. These ideas were used by Kostin to derive a 
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dissipative nonlinear Schrodinger equation [13], and Al- 
brecht [ll| worked out the general principles limiting the 
acceptable form of the nonlinear terms. 

The problem of dissipation in quantum mechanics is 
closely connected to the use of mixed quantum-classical 
treatments [l^.fl3| whereby only a few degrees of freedom 
of a complex system are treated quantum-mechanically 
while the remaining ones are described classically. Al- 
though Egorov and coworkers [13] point out that not 
always an accurate description is attained and so this 
approach must be used with care, it is recognized as a 
method of choice when a significant reduction of com- 
plexity is desired. 

In this paper we present a derivation of a rather simple 
nonlinear Schrodinger equation that includes a dissipa- 
tive term. This equation is obtained from a consideration 
of a physical model: a quantum oscillator connected to 
an infinite classical string that acts as a means to trans- 
fer energy away from the quantum system (see also [15l]). 
The obtained Schrodinger-like equation is illustrated by 
considering two specific examples. The first one, the de- 
cay of a quasi-stationary state, is an auxiliary one and 
shows that the results obtained by using our analysis do 
not contradict the ones derived by the usual and more 
complex methods using the density-matrix formalism and 
a double set of variables. The second example, the non- 
linear resonance, is interesting on its own and demon- 
strates that the quantum-mechanical nonlinear resonance 
is physically richer than its classical analog. 

Our paper is organized in the following way. In 
Sect. II, the model system is introduced and the dissi- 
pative Schrodinger-like equation for the quantum system 
is derived. In Sect. Ill, the obtained equation is simpli- 
fied employing the rotating-wave approximation. In the 
next two Sections, some illustrations of the developed 
technique are presented: in Sect. IV we consider the de- 
cay of the quasi-stationary state, and Sect. V is devoted 
to the non-linear resonance. Finally, our conclusions are 
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formulated in Sect. VI. 



II. INTERACTION OF THE QUANTUM 
OSCILLATOR WITH A CLASSICAL CHAIN 

Let us begin by introducing our model of weakly cou- 
pled quantum and classical subsystems. The quantum 
system is modeled as a harmonic oscillator of mass m and 
spring constant k, described by the quantum-mechanical 
position and momentum operators x and p. The role of 
the classical system is played by a semi-infinite classical 
string which, for the sake of convenience, is discretized 
and represented as a chain of identical balls of mass M 
interconnected by springs of equilibrium length a and 
spring constant K. The two subsystems are linked by a 
different spring (spring constant k) connecting the quan- 
tum oscillator to the ball sitting at the end of the chain 
as depicted in Fig. [TJ The corresponding Hamiltonian 
reads 
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with x n denoting the coordinate of the n-th ball mea- 
sured from its equilibrium position, and p n the respective 
momentum. 
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FIG. 1: The layout of the model system. 

Following [16], we assume that the quantum oscillator 
is described by its own wave function ^(x, t) which solves 
the Schrodinger equation 
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with the above Hamiltonian {1} while the dynamical vari- 
ables of the balls in the chain obey the classical Hamilton 
equations 
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with the Hamiltonian operator replaced by its quantum- 
mechanical average over the state of the quantum oscil- 
lator 



(H) 



dx^*(x,t)H^{x,t). 



(4) 



In this way, quantum variables and operators do not enter 
the equations for the classical degrees of freedom and one 



obtains a consistent description. This averaging consti- 
tutes the main assumption of the quasiclassical approxi- 
mation and was used earlier to treat the coupling of clas- 
sical and quantum degrees of freedom [16]. We give some 
further justification of this crucial step below. 

Eqs. lead to the standard classical equations for the 
chain 



Mx + K(x — Xx) + kxq 
Mx n - K(x n _ x - 2x n + x n+ i) 



k(x), (5a) 
0, n > 1. (5b) 

The second of these equations (f5b|) is readily solved by 
the Fourier transform 
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producing the standard dispersion law 
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From here on we restrict our consideration to the long- 
wavelength approximation (qa « 1) thus neglecting the 
discretization and effectively returning to a continuous- 
string model for the classical subsystem. In this limit, 
we have a linear dispersion law 
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w ith th e wave-propagation velocity v = Qa and Q = 

Inserting the Fourier transform ((6]) into the boundary 
condition (|5a|) and using the long- wavelength approxima- 
tion we are able to simplify its left-hand side to 
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The resulting equation of motion for the first ball coor- 
dinate Xq 
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(10) 



may be interpreted as the equation of motion for a 
dissipative classical mode. Eq. lflO|) together with the 
Schrodinger equation ^ for the quantum variables con- 
stitute the complete equation set describing the behavior 
of the dissipative quantum system. 

Now we turn to the examination of the validity of the 
employed quasiclassical approximation. A smooth way 
to proceed is to introduce the dimensionless variables by 
means of scaling 
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Here I and Iq are, respectively, the characteristic length 
scales of the quantum and classical subsystems, and Uq 
stands for the frequency of the quantum oscillator. The 
scaling enables us to rewrite the obtained set of equations 

as 

&^(x,t)=H*(x,t), (12a) 

H =\f + \& ~ A ^o, (12b) 

± + A£a:o = A(a;), (12c) 

with 

K lo 

a "^V^m' (13) 

Note that we excluded the terms which do not depend 
on the operators x and p from the quantum Hamilto- 
nian i|12b| as they can always be absorbed into the wave- 
function phase. The obtained equations describe the in- 
teraction of a quantum-mechanical system with a single 
dissipative classical mode xq. In fact, the derivation does 
not depend on the assumption that the quantum system 
is a harmonic oscillator. Thus, the equations are general 
and can be applied to any quantum system with dissipa- 
tion. For instance, considering an anharmonic oscillator 
driven by an external time-periodic force (we treat this 
case as an illustration) the above Hamiltonian is supple- 
mented with the following additional terms 

AH = -xf cos(urt) + ax 4 , (14) 

with / expressing the force amplitude and a being the 
anharmonicity coefficient. 

In order to have an appreciable (that is, not negligi- 
bly small) interaction between the quantum and classical 
subsystems one has to assume the characteristic frequen- 
cies loo and Q to be of the same order of magnitude. Un- 
der this assumption the quantum oscillator will be able 
to emit "phonons" into the string. Furthermore, the cou- 
pling constant n representing the interaction of quantum 
oscillator with the chain should not exceed the constant 
k describing the oscillator potential itself. Otherwise, 
the oscillator cannot be considered as a separate entity. 
In view of the above constraints, the dimensionless cou- 
pling constants A and £ depend essentially only on the 
adiabatic parameter m/M. 

Thus, in the case of a light quantum oscillator inter- 
acting with a chain of heavy balls the above adiabatic 
parameter is small and so are the coupling constants. 
The smallness of the coupling constant ensures a weak 
perturbation of the chain. The chain under consider- 
ation may be treated as a collection of noninteracting 
harmonic modes. Meanwhile, the quantum-mechanical 
description of a harmonic mode gives the same result as 
the classical one. Therefore, one may conclude that the 
smallness of the adiabatic parameter m/M guarantees 
the validity of the applied quasiclassical approximation 



for the description of the interaction between the quan- 
tum and classical subsystems. A similar argument was 
also presented in the development of the mixed quantum- 
classical dynamics [l2T |. 

Alternatively, one may argue that due the smallness of 
the adiabatic parameter the characteristic length of the 
wave functions of heavy balls is much smaller than the 
characteristic length I of the wave function of the oscil- 
lator. Consequently, the balls may be considered to be 
following their well-defined classical trajectories whereas 
the oscillator has to be described quantum-mechanically. 

III. ROTATING- WAVE APPROXIMATION 

Let us now proceed to some illustrations of the devel- 
oped formalism based on a weakly nonlinear oscillator. 
Namely, we assume that the Hamiltonian of the quantum 
subsystem consists of a sum of Eq. (|12b|) and Eq. (f!4|) . 
Besides, we restrict the treatment to the case of a weak 
nonlinearity, that is, we assume that the anharmonicity 
a, the force amplitude /, and the coupling constants A 
and £ are small quantities compared to the unity. We also 
restrict the frequency of the driving field to the immedi- 
ate vicinity of the resonant frequency of the harmonic 
oscillator which is unity in the dimensionless units. Thus 
we are able to treat the frequency deviation from the 
resonance as yet another small parameter r\ = uj — 1, 

We expand the sought wave function into the series of 
harmonic-oscillator eigenfunctions 

oo 

*(x,t) = «n(t)c^ E "VnW. (15) 

n=0 

Here, the dimensionless energies and the corresponding 
eigenfunctions read 

E n = n + ~, rf, n (x) = - L c- x2/2 H n {x), (16) 

^ y / 2 n n\^/TT 

and the symbol H n (x) stands for the n-th Hermite poly- 
nomial. Note that instead of the usual time dependence 
~ e~^ Ent of a harmonic-oscillator eigenfunction Eq. lfT5|) 
features a slightly modified exponential factor e - RujE ^ t 
that absorbs a part of the time dependence of the expan- 
sion coefficient a n (t). 

Using the above expansion of the wave function lfl5|) 
the coordinate average can be presented as 

oc 

n,m=0 (17) 

= V2Re (e fiwt a) 

where 

oo 

a = a n a n+1 y / n+ 1 (18) 

n=0 
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is the usual polarization of the quantum oscillator. 

Inserting the coordinate average lfl7|) into Eq. i|12cp . 
and solving it to the leading term in the small coupling 
constants we obtain the following expression for the co- 
ordinate of the first ball 



x 



V2X 



eft Re (e fiwt a) 



= V2Xlm (e SuJt a) 



(19) 



Next, inserting the obtained expression Ifl9|) and the 
wave function lfl~5)l into Eq. (|12a|) . projecting this equa- 
tion onto the n-th state, and neglecting all oscillating 
terms (this amounts to the usual rotating-wave approxi- 
mation) we arrive at the final set of equations for the time 
evolution of the wave function expansion coefficients 



fia„ = - {r]E n 



l \n)}a n 



-y/n{F + fi7a*)a„_i - \/n + l(F ~ 8~fa)a n+1 . (20) 

Here, for the sake of simpler notation we denoted 

F = f/2V2, 7 -A 2 /2. (21) 

The obtained set of equations 1|20|1 was solved in the 
two- and three-level approximations, and the results are 
presented in the following sections. 



IV. DECAY OF THE QUASI-STATIONARY 
STATE 

In order to illustrate the validity of the proposed 
method we consider a well-known problem of the decay 
of a quasi-stationary state. 

For this purpose, we assume that the quantum system 
has only two states and there is no anharmonicity and 
no external force (F = rj = a = 0). At the beginning 
the system is prepared in the upper level, that is, the 
boundary condition reads 



t 



ciq = 0, di = 1. 



(22) 



According to Eq. (J20J) , the behavior of this simple sys- 
tem is described by the following two equations 
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-ja*ao, with a — aQa[. (23) 



These equations are readily solved analytically, and the 
solution reads 
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(24) 



These filling factors are shown in Fig. [2] by the thick solid 
(excited level) and dashed (ground level) lines. We see 
that the quantum transition takes place during a finite 
period of time proportional to 7 -1 , the reciprocal fric- 
tion coefficient. The dotted curve shows the oscillator 
polarization 



(25) 
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FIG. 2: Decay of quasi-stationary state: |di| - thick solid 
curve, |ao| 2 - dashed curve, a - dotted curve, and xo - thin 
solid curve. 



induced during the quantum transition, and the thin solid 
line shows the classical coordinate of first ball in the chain 
Xo(t) which is calculated by means of expression 1(19)1 . As 
a matter of fact, this curve represents the coordinate of 
any ball in the chain because according to Eq. ((6)) the 
coordinate of the n-th ball can be expressed as 



x n (t) = xo(t — an/v). 



(26) 



Therefore, during the decay of the quantum state a string 
excitation of the shape shown by thin solid line is emit- 
ted and travels away from the quantum oscillator with a 
constant velocity v. It can be interpreted as the classical 
analog of an emitted phonon. Naturally, the total energy 
carried by the string excitation is equal to the energy dif- 
ference between the two quantum levels. Consequently, 
we may conclude that the quantum subsystem provides 
the energy quantization in the classical one. 

It is worth mentioning that taking into account the 



norm conservation (|ao| 



01 



1) Eqs. (23) can be 



transformed into another set of equations governing the 
behavior of three real variables: the inversion / = |ai| 2 — 
I ao 1 2 and the (complex) polarization a. The linearized 
version of these equations in the vicinity of the stationary 
state n = (J = — 1) reads 



I : 



-2 7 J, a 



-7a. 



(27) 



It is remarkable that these equations coincide with those 
considered in the density matrix formalism [T3|, more- 
over, from our Eq. ((27)) we are able to extract the relax- 
ation times: the longitudinal relaxation time is T\ = l/2j 
and the perpendicular one is T2 = 1/7. We see that these 
relaxation times obey the relation T2 = 2T\ inherent to 
the pure state when no ensemble average is carried out 
and this fact confirms the adequacy of our attempt to de- 
scribe a dissipative quantum system in the Schrodinger- 
equation formalism. 



5 



Two more comments are in order. First, according to 
Eqs. (|22|) we used the boundary condition at t = — oo. 
This is related to the classical description of the chain ex- 
citations when there is no spontaneous emission. When 
solving the equations numerically one has to add some 
initial vibration field fluctuation in order to force the 
quantum transition. 

Another point is related to the Fourier transform of 
the coordinate of the first ball 

*oM = r due^xoit) - (28) 

J-oo cosh(?77r/27) 

which gives us the shape of the emission line. We see 
that close to the resonant frequency, when \tj\ <§C 1, the 
line shape is Lorentzian coinciding with the result fol- 
lowing from a simple calculation based on the Fermi's 
golden rule. However, further away from the resonance, 
the line shape deviates from the Lorentzian and features 
exponential tails. 

V. RESONANT POWER ABSORPTION 

The following problem is that of the resonant power 
absorbtion which we calculated solving the equation set 
([20)1 using two-level 

do = fi|a + 8(F- 670)01, (29a) 

di = fiyfli +&(F + 8ja*)a , (29b) 
a = aoa*. (29c) 

and three-level 

d = fi (tj/2 - a) a + fi(F - 670)01, (30a) 
di = fi (377/2 — 5a) ai 

+8(F + 8 7 a*)a + V2fi(F-fi 7 a)a 2 , (30b) 
d 2 = fi (577/2 - 13a) 01 + V213(F + fi 7 a*)ai, (30c) 
a = ao^i + v2aia2. (30d) 

approximations. Note that we excluded a nonlinear term 
from the equations of the two-level approximation be- 
cause in this case it leads only to an inessential shift of 
the resonant frequency. The power absorption was cal- 
culated by averaging the instantaneous power over the 
period of the external force, that is, 

P = -2FIma. (31) 

The power absorption in the two-level approximation 
as a function of the deviation 77 is shown in Fig. [3] for 
various force amplitudes. Note that Eqs. ((29j) and f30|) 
depend linearly on their parameters, therefore, by us- 
ing an appropriate scaling of the time we may reduce 
the number of parameters by one. Therefore, we choose 
to express the quantities F, P, 77 and a in 7 units. 



We observe the following typical behavior: For a weak 
force (Fq < 0.5 7 ) the shape of the resonance resembles 
a Lorentzian, and its peak value is increasing with in- 
creasing force. When the force exceeds the critical value 
Fq = 0.5 7 the peak stops growing and is flattened at the 
top due to the saturation at the critical power absorption 
value Pq = 0.5 7 indicated by the thin dotted line. 
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FIG. 3: Power absorption in the two-level approximation. 
The numbers on curves indicate the values of the external 
force amplitude F/'y. 



In the next two Figures the results for power ab- 
sorption calculated in the three-level approximation are 
shown. In Fig. [4] we plot the power absorption in a linear 
oscillator (a = 0). 
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FIG. 4: Power absorption in the three-level approximation in 
the case of the linear oscillator. The numbers on the curves 
indicate the values of the external force amplitude F/'y. 



The general shape of the absorption curves closely 
resembles those obtained in the two-level approxima- 
tion, except that now the saturation sets in at a larger 
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force value Fq — 7, and the saturated power absorption 
Po = 1.57 is larger as well. 

A more essential difference is seen in the uppermost 
curve corresponding to the force whose amplitude (F = 
57) strongly exceeds the critical value. Here we see a 
narrow gap in the absorption which appears close to the 
resonant frequency. It resembles an analogous gap that 
appears in the power absorption calculated by means of 
the density matrix equation technique [l8[ and is caused 
by the interaction of coherent light with an inhomoge- 
neously broadened resonance line. 




of the quantum anharmonic oscillator in the absorption 
gap more carefully. In Fig. [6] the right-hand side of this 
absorption gap for the parameter values F = O.57 and 
a = 0.257 is shown in more detail. The simplest way to 
distinguish various types of dynamic behavior is to plot 
the phase portrait of the system. As there are three com- 
plex variables constrained by a single normalization con- 
dition the phase space is five-dimensional. Therefore, we 
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FIG. 6: Power absorption gap for the parameter values F 
O.57 and q = O.257. 



FIG. 5: Power absorption in the three-level approximation in 
the case of an anharmonic oscillator obtained for the param- 
eter F = O.57. The curves are numbered by the values of the 
anharmonicity coefficient a/7. 



restrict ourselves to a projection of the phase space onto 
a two-dimensional plane (Reao,Reai) shown in Fig. [7l 
The panel numbers correspond to the frequencies indi- 
cated by dots in Fig.[6l We see that on the regular flank of 



In Fig. [5] the results for the nonlinear oscillator are 
presented for a set of values of the anharmonicity coeffi- 
cient and the force amplitude corresponding to the criti- 
cal value of the two-level approximation, i. e. F = O.57. 
We see that for the linear-oscillator case (a = 0) the 
absorption demonstrates the expected Lorentzian behav- 
ior. When the anharmonicity is incremented the reso- 
nant peak moves to the right, becomes asymmetric and 
the saturation manifests itself as a widening absorption 
gap close to the shifted resonance frequency. 

The absorption gap resembles the gap in nonlinear cy- 
clotron resonance which was described using the bal- 
ance equations. In that work the gap was related to 
the dynamical chaos appearing in the classical nonlin- 
ear equations of motion. Quantum chaos is a subject 
of great interest as well (see [IS] for references), and is 
largely concerned with the quantum signatures of classi- 
cally chaotic systems. Such signatures are often sought 
by comparing the wave-function patterns corresponding 
to various eigenstates with the classical trajectories [2l| . 

It follows from our quasiclassical consideration, 
Eqs. l(2fj|) . that the behavior of the quantum system can 
also be visualized and interpreted in terms of phase-space 
trajectories obtained by treating the wave function ex- 
pansion coefficients a n as generalized coordinates. In or- 
der to illustrate this possibility, we looked at the behavior 
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FIG. 7: Two-dimensional projection of the phase space. The 
numbers on the panels correspond to the frequencies indicated 
by dots in Fig. [6] 



the resonance curve (point 1) the oscillator demonstrates 
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a rather simple behavior — the limit cycle. As a mat- 
ter of fact, in this case all three coefficients a n oscillate 
with the same frequency, namely, a n (t) = g n cxp(6(t). 
Therefore, in this case the power absorption problem can 
be simplified essentially and transformed into a nonlin- 
ear eigenvalue problem defining the frequency £, and the 
stationary coefficients g n are simply related to the filling 
factors of the harmonic oscillator levels. 

The behavior of the anharmonic oscillator in the ab- 
sorption gap (points 2, 3, 4) is quite different. Here, the 
coefficients a n oscillate with multiple frequencies, more- 
over, the oscillator can demonstrate both periodic be- 
havior (point 2) and non-periodic behavior (points 3, 4). 
This non-periodicity is in fact responsible for the gap in 
the power absorption as it causes the loss of coherence 
between the driving force and the system response. This 
is the so-called quasiperiodic behavior but the appear- 
ance of a strange attractor is also expectable at larger 
values of the anharmonicity coefficient a. 

VI. CONCLUSIONS 

In conclusion, we propose a Schrodinger-like equation 
for the description of a dissipative quantum-mechanical 



system. This equation is derived on the basis of a qua- 
siclassical treatment of coupled quantum and classical 
subsystems and justified by the difference of masses of 
the two subsystems which leads to the smallness of adi- 
abatic parameters. The considered decay of a quasi- 
stationary state demonstrates that the proposed descrip- 
tion gives identical results to those obtained from the 
density-matrix approach when only homogeneous broad- 
ening is taken into account. The model calculation of the 
nonlinear absorption indicates that many classical phe- 
nomena — such as the asymmetry of the resonance shape 
and absorption dips — brought about by the loss of co- 
herence between the driving force and the response may 
be featured by a quantum system as well. We believe 
that the proposed Schrodinger-equation formalism with 
dissipation may find application to quantum nanostruc- 
tures and quantum chaos. 
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